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LONG-TERM  GOALS 

The  long-term  goals  of  this  research  are  to  understand  the  predictability  of  acoustic  variability  in  both 
deep  and  shallow  water  ocean  environments. 

OBJECTIVES 

The  objectives  of  this  work  are  the  development  of  accurate,  and  computationally  efficient,  reduced- 
physics  acoustic  propagation  models  for  the  prediction  of  the  statistics  of  ocean  acoustic  signals.  Three 
potentially  useful  models  are  being  considered;  1)  A  coupled  mode  theory  for  acoustic  propagation 
through  random  fields  of  sound  speed  perturbations,  like  those  from  the  Garrett-Munk  (GM)  internal 
wave  spectrum,  or  modification  there-to  for  shallow  water,  2)  A  coupled  mode  theory  for  acoustic 
propagation  through  shallow  water,  non-linear  internal  waves,  and  3)  Miink  and  Zachariasen  (Rytov) 
theory  for  the  prediction  of  spectra  of  phase  and  log-amplitude  of  low  frequency,  relatively  short 
range  propagation  along  a  geometric  ray  path. 

APPROACH 

The  approach  to  this  work  is  to  first  test  the  aforementioned  models  using  numerical  experiments,  to 
establish  the  regimes  of  validity  (i.e.  range,  frequency,  acoustic  path,  etc)  with  known,  ocean  sound 
speed  fluctuation  models.  Once  the  acoustic  scattering  theory  has  been  tested  in  this  way,  the  models 
can  be  used  to  interpret  data,  where  the  environmental  conditions  have  some  degree  of  uncertainty. 
Experiments  which  will  be  used  for  model  testing  are,  1)  the  recent  Long  Range  Acoustic  Propagation 
Experiment  (LOAPEX),  and  SPICE04  data,  2)  the  North  Pacific  Acoustic  Laboratory  (NPAL)  data,  3) 
the  Acoustic  Thermometry  of  Ocean  Climate  (ATOC)  data,  4)  the  South  China  Sea  data,  5)  the 
Shelfbreak  PRIMER  data,  6)  the  SLICE89  data,  and  7)  the  SW06  data. 

WORK  COMPLETED 

Work  completed  in  the  previous  year  has  focused,  by  enlarge,  on  the  numerical  testing  of  the  three 
aforementioned  models.  Regarding  the  coupled  mode  model  for  propagation  through  random  sound 
speed  fields,  we  have  tested  the  model  in  deep  water  settings  at  a  single  frequency,  and  we  have 
worked  on  the  incorporation  of  broadband  terms  into  the  theory.  The  broadband  model  still  needs  to  be 
validated  by  numerical  experiment.  The  coupled  mode  model  for  shallow  water,  non-linear  internal 
waves  was  formulated  and  tested  by  numerical  experiment  over  the  year,  and  we  are  ready  to  work  on 
comparisons  to  observations.  The  Munk  and  Zachariasen  model  has  had  some  numerical  testing  and 
has  been  compared  to  observations  from  the  ATOC  87-km  range,  75  Hz  frequency  Acoustic 
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engineering  Test  (AET).  More  extensive  numerieal  testing  of  the  Munk  and  Zaehariasen  model  are 
underway  right  now,  examining  the  regimes  of  applieability  of  the  model  as  a  function  of  frequency, 
range,  and  acoustic  path. 

RESULTS 

A,  Coupled  Modes:  Sound  Transmission  through  Stochastic  Ocean  Internal  Waves 

Audrey  Morozov  (of  Woods  Hole  Oceanographic  Institution)  and  I  have  worked  to  develop  a  coupled 
mode  theory  to  better  understand  acoustic  intensity  statistics.  The  approach  is  largely  based  on  the 
work  of  Dozier  and  Tappert  (1978),  and  of  Van  Kampen  (1992),  but  our  approach  modifies  their 
techniques  enabling  the  computation  of  the  important  cross  mode  coherences.  Broadband  effects  are 
also  treated  by  computing  the  cross  mode  coherences  across  frequency.  The  primary  mode  coupling 
mechanism  in  the  theory  is  a  Bragg-like  resonance  between  the  difference  modal  wave-number  and  the 
internal  wave  wave-number.  This  work  closely  parallels  efforts  by  Alex  Veronovich  and  Frank 
Henyey.  We  also  intend  to  extend  the  theory  to  deal  with  shallow  water  environments  in  which  there  is 
significant  modal  attenuation. 


Sound  source  depth  2  km,  250Hz,  horizon  1 .3  km 


Figure  1:  Monte  Carlo  Simulation  (Red)  and  model  prediction  (Blue)  of  mean  sound 
intensity  for  250  Hz  sound  propagation  through  random  fields  of  GM  internal  waves. 

In  this  example  the  source  depth  is  2-km,  and  the  receiver  depth  is  1.3  km.  The 
unperturbed  sound  intensity  is  plotted  in  green  for  reference.  The  model  prediction  is 
good  to  within  a  few  dB,  with  biggest  differences  at  low  intensity.  The  unperturbed 
calculation  for  the  mean  intensity  is  clearly  not  appropriate. 


Initial  comparisons  between  our  model  and  Monte  Carlo  parabolie  equation  simulations  through 
Garrett-Munk  internal  waves  at  250  and  100  Hz  show  the  theory  predicts  the  multi-megameter  range 
evolution  of  the  mean  intensity  to  within  a  few  dB.  The  theoretieal  caleulation  takes  only  a  few 
minutes  to  eompute  while  the  Monte  Carlo  ealeulation  takes  several  days.  An  example  of  the  mean 
intensity  computed  from  the  theory  and  numerieal  simulation  is  shown  in  figure  1 .  The  agreement 
between  the  model  and  the  Monte  Carlo  simulations  is  exeellent. 

This  year  we  have  worked  hard  to  include  broadband  effeets  into  the  model  by  eomputing  eross  modal 
coherenees  aeross  different  frequencies  (Figure  1  only  shows  the  single  frequeney  result).  The 
broadband  theory  has  been  worked  out,  but  due  to  memory  limitations  on  our  present  workstations,  we 
have  not  been  able  to  eompute  direet  eomparisons  with  the  broadband  Monte  Carlo  calculations.  We 
expeet  to  have  new  equipment  soon  to  finish  the  ealeulation. 

B,  Coupled  Modes:  Shallow  Water,  Non-linear  Internal  Waves 

Development  of  a  model  for  shallow  water  propagation  through  non-linear  internal  waves  has  been 
very  sueeessful  this  year.  Mueh  of  this  work  was  done  as  part  of  the  NFS  masters  thesis  earned  out  by 
Lt  Aaron  Young  (AUS).  Using  methods  from  quantum  mechanieal  time-dependent  perturbation 
theory,  the  following  expression  was  obtained  for  the  ehange  in  mode  energy  after  acoustic 
propagation  through  a  single  internal  solitary  wave  (ISW)  of  width  A,  position  ro,  and  maximum 
fraetional  sound  speed  po, 

I  A(R)  t=l  Am  r  +2A;r''V.Zl,/f,(0)/(;(0)C.,.  exp(-/iAV4)sm(Z,.,r,) 

where  |An(0)|  ,  |An(R)|  is  the  energy  in  mode  n  before/after  propagation  through  the  ISW,  Cmn  is  a 
coupling  matrix  involving  the  acoustie  mode  funetions  m  and  n  and  the  depth  structure  of  the  ISW,  and 
Imn  is  the  beat  wavenumber  between  modes  m  and  n.  The  physieal  model  here  is  of  a  Bragg  type 
resonance  between  the  ISW  horizontal  wavenumber  component,  that  matehes  the  beat  wavenumber 
Imn-  It  is  important  to  note  that  although  a  perturbation  theory  approaeh  was  used  here  the  model  does 
eonserved  energy.  Further,  the  model  ean  also  be  used  to  prediet  mode  energy  ehanges  through  ISW 
packets  or  solibores,  and  it  ean  be  used  to  prediet  other  aeoustie  observables  like  coherenee  or 
seintillation  (more  later).  Figure  2  shows  how  the  model  predietions  compare  to  parabolie  equation 
simulations  as  the  ISW  width  is  ehanged  from  20-m  to  1000-m.  The  model  is  elearly  better  at  lower 
frequeney,  but  it  is  not  eompletely  inaccurate  at  the  highest  frequeney  tested  (300-Hz).  As  expeeted, 
coupling  is  seen  to  be  small  for  very  narrow  and  very  wide  ISWs,  and  maximum  coupling  is  obtained 
for  an  ISW  width  roughly  equal  to  the  beat  wavenumber  Imn. 
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Figure  2:  Parabolic  equation  (Blue)  and  model  estimates  (Red)  of  mode  energies  after 
propagation  through  an  ISW  with  a  maximum  sound  speed  perturbation  of  15  m/s  for 
various  ISW  widths  A  from  20-1000  (m).  The  upper,  middle,  and  lower  rows  are  for 
acoustic  frequencies  of  75,  200,  and  300  Hz.  The  straight  dash  curves  show  the  initial 
mode  energy  at  the  source.  The  mode  numbers  chosen  in  the  display  are  the  four 
highest  energy  modes  excited  by  the  source. 


Also  of  interest  are  the  model  predietions  of  mode  energy  as  a  function  of  ISW  position  ro.  Figure  3 
shows  the  variation  of  mode  energy  for  ISW  positions  between  1000  and  5000-m  for  the  150  Hz 
frequency,  and  an  ISW  width  of  lOO-m.  The  agreement  is  seen  to  be  excellent.  Figure  2  (right  panel) 
also  shows  results  for  the  case  of  coupling  through  a  packet  of  three  ISWs  of  strength  15  m/s,  13.5 
m/s,  and  12  m/s;  again  the  agreement  is  seen  to  be  excellent.  Note  that  the  pattern  of  mode  energy 
variation  as  a  function  of  ISW  position  is  not  regular  or  repeating,  but  is  modulated  by  the  non- 
commensurate  modal  wavenumber  difference  Imn. 


Work  is  already  underway  to  evaluate  second  order  terms,  which  will  improve  the  model  at  higher 
frequency.  In  addition  the  model  is  being  adapted  to  predict  other  acoustical  observables  like 
coherence  and  scintillation. 


Figure  3:  Parabolic  equation  simulations  (solid)  and  model  predictions  (dash)  of  mode 
energies  after  propagation  through  ISW’s  of  width  100-m,  at  a  frequency  of  150-Hz, 
for  ISW  positions  from  1000-5000  m.  The  seven  curves  correspond  to  the  seven  most 
energetic  modes.  In  the  left  panel  a  single  ISW  is  modeled  having  a  maximum  sound 
speed  perturbation  of  15  m/s,  while  the  right  panel  models  a  three  ISW  packet  with 
maximum  sound  speed  perturbations  of  15, 13.5,  and  12  m/s  for  each  of  the  three 
ISWs.  The  agreement  between  model  and  numerical  experiment  is  excellent. 

The  variation  of  mode  energy  as  a  function  of  ISW  position  in  this  2-D  model  will  have  important 
consequences  for  situations  in  which  the  ISW  propagation  is  at  some  small  angle  relative  to  the 
acoustic  propagation  path.  In  this  case  fully  3-d  effects  (i.e.  horizontal  refraction)  are  small,  and  the 
present  model  can  be  used  in  an  Nx2D  capacity.  Two  effects  occur  when  there  is  a  tilted  ISW;  First  the 
effective  distance,  ro,  to  the  ISW  changes  as  a  function  of  azimuthal  launch  angle  from  the  source,  and 
second  the  ISW  appears  to  have  a  large  width.  A,  as  the  mode  intersects  the  ISW  at  a  more  glancing 
angle.  Figure  4  shows  just  such  a  calculation  for  a  tilted  ISW.  The  striation  pattern  in  the  mode  energy 
after  propagation  through  the  ISW  is  primarily  caused  by  the  variation  of  ro,  along  the  ISW,  and  not 
the  effective  change  in  the  ISW  width.  Work  on  this  3-d  problem  was  part  of  an  NFS  masters  thesis 
composed  by  Lt  Doug  Roush.  Further  numerical  testing  of  this  Nx2D  application  of  the  model  will 
need  to  be  carried  out  to  establish  what  range  of  ISW  angles,  packet  parameters,  and  acoustic 
frequencies  in  which  out  of-plane  scattering  is  small. 
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Figure  4:  Model  prediction  for  the  ratio  of  unperturbed  to  perturbed  mode 
1  energy  for  300-Hz  acoustic  propagation  through  a  tilted  solution  (black  line).  The 
point  source  is  located  at  the  origin.  The  striation  pattern  past 
the  ISW  is  due  primarily  to  the  variation  of  ISW  position  along  the 
ISW front,  rather  than  the  apparent  broadening  of  the  ISW from 

the  glancing  angle. 


C,  Weak  Fluctuation  Theory 

A  PhD.  student  in  the  WHOI/MIT  Joint  program,  Mr  Jinshan  Xu,  and  I  are  applying  the  weak 
fluctuation  theory  of  Munk  and  Zachariasen  (1976)  to  the  prediction  of  frequency  and  vertical 
wavenumber  spectra  of  phase  and  log-amplitude  for  acoustic  propagation  through  random  sound  speed 
fields  like  those  described  by  the  GM  spectrum.  Here  the  spectra  are  written  as  integrals  along  a 
geometric  acoustics  ray  path,  and  the  acoustic  fluctuations  are  only  influenced  by  internal  wave 
wavenumbers  that  are  perpendicular  to  the  ray  path  (assumed  locally  straight).  Accordingly  the 
frequency/wavenumber,  (co,kz)  spectrum  of  phase  (j),  and  log-amplitude  y,  at  range  R  along  ray  path  T 
is. 


j  ds^{k^  (co,  );  ) 


l±cos 


KR}^ 

2n 


2 

where  ko  is  a  reference  acoustical  wavenumber,  Rf  is  the  first  vertical  Fresnel  zone,  N  is  the  buoyancy 
frequency  profile,  and  O  is  the  spectrum  of  internal  waves  evaluated  at  the  perpendicular  wave  number 
given  by 


2  2  2  1/2 

Note  here  that  col  =(f  +N  (zray)  tan  (Gray))  ,  (f  is  the  loeal  Coriolis  parameter),  is  a  low  frequeney  out 
off  that  depends  on  the  ray  angle;  thus  steep  rays  will  not  be  able  to  interaot  with  low  frequeney 
internal  waves.  Further  there  is  a  high  frequeney  out  off  at  the  loeal  buoyanoy  frequeney  N,  thus  the 
speotrum  has  two  Heavyside  funotions  H  to  impose  the  low  and  high  frequeney  outoffs. 


This  analysis  was  motivated  by  the  1994  Aooustio  Engineering  Test  (AET)  data  whioh  demonstrated 
very  weak  aooustio  fluotuations  of  75  Hz  sound  transmissions  to  87-km  range.  Eigure  5  shows  the 
oomparison  between  Monte  Carlo  parabolio  equation  simulation  through  random  fields  of  Garrett- 
Munk  internal  waves,  and  the  theory  for  parameters  olose  to  that  of  the  AET.  The  model  is  seen  to 
olosely  prediot  the  speotra  even  when  a  broadband  signal  is  used  (The  Munk  and  Zaohariasen  theory 
assumes  a  narrowband).  Presently  anNPS  master’s  student,  Et  Serdar  Tombul  (TK)  is  working  on 
extending  the  model  oomparisons  to  frequenoies  of  75,  150,  and  300  Hz,  4  ray  path  geometries  from 
axial  to  surfaee  refieeting,  and  upto  ranges  of  200  km. 
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Figure  5:  Frequency  spectra  of phase  and  log-amplitude  from  broadband  (red  dot), 
and  narrowband  (dash)  parabolic  equation  simulations  of  sound  transmission  to  87- 
km  range  through  random  fields  of  Garrett-Munk  internal  waves.  Black  curves  show 
the  Munk  and  Zachariasen  (Rytov)  predictions.  Two  arrivals  are  modeled:  one  has  a 
wavefront  identifier  of  -3,  and  the  other  has  a  value  of  +4. 

Next  we  compared  the  model  to  the  log-amplitude  frequency  spectra  from  the  AET.  The  comparison  is 
limited  to  only  relatively  high  and  low  frequencies  because  of  the  sparse  time  sampling  of  the  AET, 
however  Eigure  6  shows  that  the  agreement  between  model  and  observations  is  satisfactory.  This 
comparison  explains  the  rather  curious  observational  result  that  the  wavefront  ID,  -3,  which  traverses 
the  sound  channel  at  a  higher  angle  has  much  less  low  frequency  variability  than  the  ID  which  has 
smaller  angles  and  does  not  get  as  close  to  the  surface.  This  attenuation  of  low  frequency  variability 
for  steep  rays  is  the  consequence  of  the  resonance  condition  between  the  local  ray  tilt  and  the  internal 
waves  with  wave  numbers  perpendicular  to  this  tilting  ray.  We  believe  that  this  result  is  the  first 
observational  evidence  for  this  resonance  condition. 


Figure  6:  Observed  and  modeled  frequency  spectra  of  log-amplitude  from  the  87-km 
range,  75  Hz,  AET.  The  observed  spectra  for  wavefront  ID’s  -3  (blue  dots),  and  +4 
(black  dots)  should  be  compared  to  the  predictions,  red  for  ID  -3,  and  black  for  ID  +4. 

The  gap  in  the  observed  spectra  is  due  to  irregular  sampling. 

IMPACT/APPLICATIONS 

There  are  several  implieations  of  this  work  to  the  understanding  of  acoustic  predictability.  A  short  list 
of  the  major  issues/impacts  are  given  below. 

1)  Our  broadband  approach  to  coupled  mode  theory  for  sound  transmission  through  random  fields 
of  internal  waves  may  provide  the  first  theoretical  framework  to  understand  the  in-filling  of 
shadow  zones  from  acoustic  scattering.  Examples  include  the  vertical  broadening  of  the  wavefront 
finale,  the  vertical  extension  of  wavefront  branches  observed  at  Navy  SOSUS  stations,  and  the 
energy  scattering  into  the  sound  channel  axial  region  from  off-axis  sources.  The  theory,  so 
modified  to  include  modal  attenuation,  may  also  have  application  in  shallow  water  environments. 

2)  Progress  on  a  coupled  mode  theory  for  low  frequency  long  range  propagation  may  be  able  to 
fill  the  gap  in  understanding  left  when  the  path  integral  theory  of  acoustic  fluctuations  was  shown 
to  be  seriously  flawed  past  propagation  ranges  of  a  few  hundred  kilometers. 

3)  Our  2-d  coupled  mode  model  for  sound  transmission  through  single  or  multiple  ISWs  is  one  of 
the  rare  examples  of  a  reduced  physics  model  for  shallow  water  acoustic  variability.  The  model 
may  be  applied  to  other  acoustic  observables  besides  mode  energy,  like  coherence  and  scintillation. 
The  model  could  also  be  used  in  a  stochastic  form  given  stochastic  models  for  all  the  ISW 


parameters  like  ISW  width,  phase  speed,  amplitude,  and  ISW  number.  Finally,  as  demonstrated 
here,  the  model  can  be  used  in  an  Nx2D  mode  to  treat  3-D  effects,  and  intensity  striation  patterns. 

4)  The  Munk  and  Zachariasen  theory  will  allow  us  to  quantify  the  space-time  scales  of  acoustic 
variability  for  low  frequency,  sound  propagation  over  a  few  convergence  zones.  This  information 
leads  directly  to  models  of  coherence  that  can  be  exploited  to  construct  "optimal”  array  and  signal 
processing  algorithms. 
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